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(54) Abstract Title 

Correcting defective pixels In an Image 



(57) A source imago formed from a plurality of pixels 
each having a respective value may have defective pixels 
(X}. A defective pixel (X) is corrected based on curvature 
Information computed from pi*el values of pixels located 
near the defective pixel (X). Curvature is defined as being 
calculated by performing arithm&tic functions on values of 
pixels in rows above and balow a row containing the 
defective pixel (X) or on values of pixels In columns 
preceding and succeeding the column containing the 
defective pixel (X). A second Independent claim is 
Included which claim* the correction of defective pixel 
values by determining a median pixel value from pixels on 
the same colour plane located near the defective pixel (X). 
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CORRECTING DEFECT PIXELS IN A DIGITAL IMAGE 

TECHNICAL FIELD ? 

This invention relates to systems and methods for correcting defect pixels 
in a digital image (for example, a digital still image, a digital video image, or 
5 both), 

BACKGROUND 

Digital images are formed from an array of pixels (picture elements), each 
pixel having an assigned color and intensity. A digital image may be captured by ' 
a charged-coupled device (CCD) with a planar array of photodetectors. Each 

10 photodetector is located at a respective pixel position. Defective photodetector* 
produce imperfections in the corresponding pixel locations in the resulting digital 
image. Although CCD manufacturers have reduced the photodetector defect rate 
significantly, elimination of all defects typically is not practical Consequently, 
manufacturers of digital imaging devices generally include systems designed to 

is correct imperfections caused by photodetector defects. Various correction 

systems have been proposed. Some correction systems include image processors 
designed to interpolate the correct values of defect pixels based upon the values 
of non-defective neighboring pixels, Such correction systems may be classified as 
directional interpolation systems, or non-directional averaging (or erosion) 

20 systems. For example, in one proposed directional interpolation system, a defect 
pixel value is obtained by linear interpolation of a tine of pixels that extends 
through the defect pixel and is oriented in a direction determined from gradient 
information derived from the pixels neighboring the defect pixel. Typical non- 
directional averaging systems tend to be easier to implement; however, they also 

25 tend to overly smooth non-edge image areas and significantly distort edges in the 
corrected image. 

SUMMARY 

The invention features improved systems and methods of correcting one or 
more defect pixels in a source image formed from a plurality of pixels each having 
30 a respective value. 



i 



In one aspect, the invention features a scheme !n which a. defect pixel is 
corrected based upon curvature information computed from pixel values located 
near the defect pixel 

Embodiments may include one or more of the following features. 
5 The defect pixel preferably i3 corrected by interpolating values of pixels 

aligned with the defect pixel in an interpolation direction determined from the 
curvature information, The interpolation direction may be determined by 
comparing curvature values computed for two orthogonal directions. The 
interpolation direction preferably corresponds to the direction having the lower 
10 curvature value. The defect pixel preferably is assigned a value computed from a" 
linear interpolation of values of pixels aligned with the defect pixel along the 
interpolation direction. The defect pixel may be assigned a value computed from 
a linear interpolation of values of two non-defective pixels located on opposite , 
sides of the defect pixel The assigned defect pixel value may be computed from 
is the operations of adding and shifting. 

The defect pixel may be corrected based upon curvature information 
computed for two different image directions. For example, a horizontal curvature 
value may be computed based upon pixel values from two groups of horizontally 
aligned pixels located on opposite sides of the defect pixel. A vertical curvature 
20 value may be computed based upon pixel values from two groups of vertically 
aligned pixels located on opposite sides of the defect pixel Alternatively, 
curvature values may be computed based upon values of pixels surrounding the 
defect pixel. The defect pixel may be corrected based upon a comparison of the 
computed curvature values, For example, the defect pixel may be corrected by 
25 interpolating values of pixels aligned vertically with the defect pixel in response to 
a determination that the horizontal curvature exceeds the vertical curvature. 
Alternatively, the defect pixel may be corrected by interpolating values of pixels 
aligned horizontally with the defect pixel in response to a determination that the 
vertical curvature exceeds the horizontal curvature. 
30 The source image may be obtained from an image sensor having a plurality 

of light detectors each generating a pixel value of a respective color, wherein the 
colors are arranged in multiple color planes of non^overlapping pixel locations. 



The defect pixel may be corrected based upon curvature information computed 
from values of pixels in the same color plane. ' ^ 

In another aspect, the invention features a scheme in which a median pixel 
value is determined from values of pixels an the same color plane located near a 
s defect pixel, and the defect pixel is corrected based upon the median pixel value. 
Embodiments may include one or more of the following features. 
The median pixel value preferably is determined from values of pixels on 
the same color plane adjacent to the defect pixel A difference value 
representative of the difference between the defect pixel value and the computed 
10 median pixel value may be computed, and the defect pixel may be corrected by 
replacing the defect pixel value with the median value in response to a 
determination that the difference value exceeds a threshold value. 
Among the advantages of tfos invention are the following. 
The curvature-based directional interpolation invention provides a robust 
is and efficient scheme for recovering corrupted pixel data in both low frequency 
(smooth) areas and high frequency (edge and corner) areas of the source image. 
The median filter invention, on the other hand, provides an efficient scheme for 
recovering corrupted pixel data on the fly, especially when the locations of defects 
are not known beforehand. The median filter Invention effectively removes 
20 outliers (i.e-, adjacent defect pixels), while substantially preserving edges and 
comers in the source image. 

Other features and advantages of the invention will become apparent from 
the following description, including the drawings and the claims. 

DESCRIPTION OF DRAWINGS 

25 FIG. 1 is a block diagram of selected components of a digital imaging 

device. 

FIG, 2 is a diagrammatic view of a pixel array with a Bayer geometry color 
pattern formed from non-overlapping color planes, 

FIGS, 3A-3D are diagrammatic views of four pixel array color planes that 
30 may be extracted from the Bayer geometry color pattern of FIG. 2. 

FIG- 4 is a flow diagram of a method of correcting a defect pixel in a digital 
source image. 



FIG. S Is block diagram of a circuit for computing a linear interpolation 
approximation. ^ 

FIG. 6 is a diagrammatic view of a median filter window superimposed on 
a pixel array. 

s FIG. 7 is a flow diagram of a method of correcting a defect pixel in a digital 

source image. 

DETAILED DESCRIPTION 

Referring to FIG. 1, a digital imaging device 10 has an imaging section 12 
and a correction filter section 14- Imaging section 12 includes a lens 16, a color 

10 filter array 18 (e.g„ a Bayer array), a CCD 20, an analog-to-digital (A/D) converter 
22, and an image buffer 24 (e.g., a random access memory). Correction filter 
section 14 includes a digital signal processor (DSP) 26 and a random access 
memory (RAM) 28. A controller 30 coordinates the operation of the components 
of imaging section 12 and correction filter section 14* Digital imaging device 10 

is may be a digital camera configured to record digital still images or a digital video 
recorder configured to record digital video images. Images captured by digital 
imaging device 10 may be recorded in a storage memory 32 (e.g., a removable 
memory card or a digital video cassette). In operation, image light is captured by 
lens 16 and directed through color filter array 18 to CCD 20, CCD 20 converts 

20 impinging filtered light into a plurality of analog image signals, and A/D converter 
22 converts the analog image signals into digital image signals. The digital image 
signals are stored in image buffer 24 as source image data for processing, After a 
sufficient amount of source image data has been stored, the source image data is 
processed by DSP 26, which corrects one or more defect pixels in the source 

25 image to produce corrected image data. The corrected image data may be 
processed further (e.g., the source image may be demosaiced), and the final 
image data may be stored in storage memory 32. 

Referring to FIG. 2, color filter array 18 is configured to transmit light of a 
single color to each pbcel location of CCD 20. For example, in a Bayer color filter 

30 array geometry, the pixel colors are arranged in multiple color planes of non- 
overlapping pixel locations. When superimposed, the color planes form a pattern 
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of chrominance colors (red and blue) interspersed among a checkerboard pattern 
of luminance colors (green), as shown in FIG. 2. ' -;j 

In the present embodiment, DSP 26 corrects defect pixels based upon 
information computed from the values of the nearest neighbor pixels; other 
embodiments, however, may incorporate into the correction process the values of 
pixels located farther away from the defect pixels. A defect pixel may be 
identified daring calibration of digital imaging device 10 using known techniques 
(e.g., filtering and testing images produced by CCD 20). The locations of defect 
pixels may be stored as a defect map in RAM 28. In the present embodiment, the 
nearest neighbors of a particular defect pixel are the eight pixels located in the 
same color plane as the defect pixel. For example, as shown in FIGS. 3A and 3B, 
in a red color plane 40 and a blue color plane 42, the eight nearest neighbors of a 
defect pixel (X, Y) are spaced apart from the defect pixel by one pixel location. . 
Referring to FIGS. 3C and 3D, two different color planes 44, 46 may be extracted 
from the Bayer geometry to identify the nearest neighbors of a defect pixel (Z). In 
color plane 44, the eight nearest neighbors are spaced apart from the defect pixel 
by one pixel location in every direction whereas, in color plane 46, four of the 
eight nearest neighbors are located adjacent to the defect pixel and the remaining 
neighbors are spaced apart from the defect pixel by one pixel location. 

As mentioned above, in one embodiment, DSP 26 corrects defect pixel 
values based upon curvature information computed from pixel values located 
near each of the defect pixels. This approach provides a robust and efficient 
scheme for recovering corrupted pixel data in both low frequency (smooth) areas 
and high frequency (edge and corner) areas of the source image. 

Referring to FIG. 4, in one embodiment, DSP 26 may correct a defect pixel 
in a source image as follows. If the curvature information is complete (i.e., none 
of the nearest neighbors is defective) (step 50), a horizontal curvature value (step 
52) and a vertical curvature value (step 54) are computed. The horizontal 
curvature (HC) and vertical curvature (VC) values may be estimated respectively 
from two groups of horizontally-aligned and vertically-aligned pixels located on 
opposite sides of the defect pixel as follows: 



Horizontal Curvature -|A-2xB + C| + |F-2xG-+H| (1) 
Vertical Curvature - |A-2*D + F1 + [ C - 2 * E + H | >^2) 



where, A, B, C, D, E, F, G and H correspond to the values of the eight nearest 
5 neighbor pixels of a defect pixel X, arranged in accordance with the following 
pixel array (i.e., A, B and C are the top neighbors, D is the left neighbor, E is the 
right neighbor, and F, G and K are the bottom neighbors). 
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10 Pixel Array I: Curvature Information is Complete 



If the horizontal curvature exceeds the vertical curvature (step 56), the defect 

pixel is corrected by interpolating values of pixels aligned vertically with the 

defect pixel (i.e., the value of X is replaced with the value of Vi * (B + C)) (step 
is 58). Otherwise, the defect pixel Is corrected by interpolating values of pixels 

aligned horizontally with the defect pixel (i,e., the value of X is replaced with the 

value of Vz*{D + E)) (step 60), 

If the curvature information is incomplete (step SO) but the gradient 

information is complete (step 62) , a horizontal gradient value (step 64) and a 
20 vertical gradient value (step 66) are computed. The horizontal gradient (HG) and 

vertical gradient (VG) values may be computed as follows; 



Horizontal Gradient - 1 1 - J | (3) 
Vertical Gradient - [ I - K 1 (4) 

25 

where I, J and K are arranged on different color planes in accordance with the 
following pixel array, 
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Pixel Array II: Incomplete Curvature Information, 
Complete Gradient Information 

s The values of other groups of pixels (e.g., A, B and D, or E, G and H), which are ' 
located at comer positions near the defect pixel, also may be used to compute the 
horizontal and vertical gradient information. If the horizontal gradient exceeds 
the vertical gradient (step 68] and the vertical interpolation information is 
complete [step 70} t the defect pixel is corrected by interpolating values of pixels 

10 aligned vertically with the defect pixel (i.e., the value of X is replaced with the 
value of Vi*(B + G)) (step 58). On the other hand, if the vertical gradient 
exceeds the horizontal gradient (step 68) and the horizontal interpolation 
information is complete (step 72), the defect pixel is corrected by interpolating 
values of pixels aligned horizontally with the defect pixel [i.e., the value of X is 

is replaced with the value of V4 • (D + E)) (step 60) , If the information required for 
the horizontal or vertical interpolation is incomplete (steps 70, 72) , the defect 
pixel value is replaced with a linear interpolation approximation in the 
interpolation direction (step 74), For example, if a horizontal interpolation is 
called for but both pixel X and pixel Y are defective, the values of pixel X and 

20 adjacent pixel Y may be recovered by linearly interpolating between pixel D and 
pixel L, as follows. 

X ~ 2/3 x D + 1/3 * L (S) 
Y = 1/3 * D + 2/3 * L (6) 

25 

where D and L are arranged on the same color plane in accordance with the 
following pixel array. 



A 




B 




C 




















D 




X 




Y 




L 
















F 




G 




Defect 







Pixel Array III: Incomplete Curvature Information, Complete 
Gradient Information, Incomplete Interpolation Information 



s The system used to compute this linear interpolation approximation may be ' 
simplified by reducing the approximation to a form that may be computed by a 
series of adding and shifting operations. For example, equations (5) and (6) 
above, may be approximated by the following equations: 

10 X * V* * D + Vi x l = VWD + VWL + D)] (7) 

Y ~ 14 x L + M x D - + Vi-(D + L)) (8) 

Equations (5) and (6) maybe approximated with greater accuracy by the 
following equations: 

is 

X ~ 5/8 * D + 3/8 * L - &*[D + W-(L + %-(D + L)J) (9) 
Y « 5/8 x l + 3/8 x D = V4-CL + V4-{D + Vi-(L + D)J) (10) 

In a preferred embodiment, the values of defect pixels X and Y are approximated 
2a by equations (9) and [10), respectively. However, equations (5) and (fi) maybe 
approximated with even greater accuracy by the following equations: 

X ~ 11/16 x p + 5/16 x l - H-(D + V4'(L + (D + x h ■ (L + D)))) (11) 
Y - 11/16 x L + 5/16 x D - Vi {L + Vi-(D + VHL + * (D + L)))) (12) 



25 



The accuracy with which equations (5) and (6) are approximated may be 
increased still further by iteratively increasing the number of adding and shifting 
operations. 
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Referring to FIG, 5, each of the approximation equations {7)-{U) may be 
computed by an approximation circuit 80. Approximation circuit 80 includes a ' Vj 
4 * 2 multiplexer (MUX) 82, an adder 84, a 1 * 2 multiplexer 86, and two registers 
88. 90. Multiplexer 82 has four inputs: two inputs are coupled to receive the two 

s pixel values (D and L) from which the defect pixel values (X and Y) are to be 
interpolated; the remaining two inputs are coupled to receive the output values 
from registers 88 and 90, respectively. A select signal 92 controls which two of 
the four multiplexer inputs are applied to the inputs of adder 84. The output of 
adder 84 - which represents the sum of the adder input values - is applied to the 

10 input of multiplexer 86. A select signal 94 controls which register 88, 9D receives 
the sum applied to the input of multiplexer 86. 

In operation, controller 30 may choreograph the operation of 
approximation circuit 80 as follows. First, through select signal 92, controller 30' 
directs multiplexer 82 to apply inputs O and L to the Inputs of adder 84, The 

15 resulting sum (L + D) is applied to register 90. Through shift signal 98, 

controller 30 shifts the contents of register 90 by one decimal place to produce the 
quantity Vi*(L + D), The quantity in register 90 is applied to the input of 
multiplexer 82. In the next cycle, controller 30 directs multiplexer 82 to apply 
inputs D and the quantity of register 90 to the inputs of adder 84, The resulting 

20 sum (D + Vz ■ (L + D]) is applied to register 88 and, through shift signal 96, 

controller 30 shifts the contents of register 88 by one decimal place to produce the 
sum (V**(D + Vi -(L + D))), which corresponds to the first approximation of the 
interpolation of defect pixel X (i.e., equation (7)). In the following cycle, the 
quantity of register 90 and input L are applied to the inputs of adder 84. The 

25 resulting sum (L + Vi - (D + L)) is applied to register 90. Controller 30 shifts the 
contents of register 90 by one decimal place to produce the quantity (Vi - (L + 
W ■ (D + L))), which corresponds to the first approximation of the interpolation of 
defect pixel Y (i.e., equation (8)). The interpolation approximations of equations 
(9) and (10] may be computed by further cycling of approximation circuit 80, For 

30 example, the sum of input D and the contents of register 90 may be applied to 
register 88, and the contents of register 88 may be shifted by one decimal place to 
produce the quantity (Vi-(D + Vi-(L + Vi-(D + L}))) (i.e., equation (9]). In 
addition, the sum of input L and the contents of register 88 may be applied to 



register 90, and the contents of register 90 may be shifted by one decimal place to 
produce the quantity (Vi*(L + Vi-(D + Vz-(L + D)))) (i.e., equation (10)). ^ 
Further cycling of circuit 80 in a similar manner would increase the accuracy with 
which the linear interpolations for X and Y may be approximated, 
s Referring back to FIG. 4, if both the curvature information (step 50) and 

the gradient information (step 62) are incomplete, the defect pixel value (X) is 
replaced with the average value computed from the values of the non*defective 
nearest neighbors (step 76); 

10 X = (P, + ? 2 + . . . + pj/N (13) 

where N is the number of non-defective nearest neighbors* and F, are the values 
of the non-defective nearest neighbors. In one embodiment, pixels are processed, 
in raster scan order. Consequently, in this embodiment, pixels above and to the 

15 left of the current defect pixel (X) are guaranteed to have valid values and 
equation (13) will produce a valid estimate of the defect pixel value. 

In another embodiment, DSP 25 corrects defect pixel values based upon a 
median pixel value computed from pixel values located near each of the defect 
pixels. This median filter approach provides an efficient scheme for recovering 

10 corrupted pixel data on the fly, especially when the locations of defects are not 
known beforehand (i.e., a defect map is not stored in memory 28). This approach 
effectively removes outliers (i.e., adjacent defect pixels), while substantially 
preserving edges and comers in the source image. 

Referring to FIGS. 6 and 7, DSP 26 may be configured to implement a 

2S median filter as follows, For each pixel value (X) (step 100), a median filter 

window 102 is defined (step 104), In this embodiment, defect pixels are corrected 
based upon information computed from the values of the nearest neighbor pixels; 
other embodiments, however, may incorporate into the correction process the 
values of pixels located farther away from the defect pixels. In this embodiment, 

30 the nearest neighbors of a particular pixel are the eight pixels (A, B t C, D, E, F 4 G 
and H) located in the same color plane as the pixel. The median value (M) of the 
values of the nearest neighbor pixels ia determined (step 106). An eight-input 
sorting (or ordering) circuit may be used to determine the median value of the 

10 



nearest neighbor pixels with the median filter window 102. If the difference 
between the pixel value (X) and the median value (M) is greater than a threshold' Vj 
value (T) (step 108) , the value of pixel X is replaced with the median value (step 
110). Otherwise, the value of pixel X is unchanged, and the next pixel is 
processed (step 100). 

The systems and methods described herein are not limited to any particular 
hardware or software configuration, but rather they may be implemented in any 
computing or processing environment. Controller 30 preferably is implemented in 
firmware, Alternatively, however, controller 30 may be implemented in a high 
level procedural or object oriented programming language, or in assembly or 
machine language; in any case, the programming language may be a compiled or 
interpreted language. 

Other embodiments are within the scope of the claims. 

For example, although horizontal and vertical directions were selected as 
the primary processing directions (i.e., the directions used to determine the 
orientation of predominate image structure and to interpolate values for defect 
pixels), other processing directions (e.g., diagonal or radial image directions), 
including two or mere directions, may be selected- 

In another embodiment, horizontal and vertical curvatures may be 
estimated by incorporating a filter array that includes three-by-three, horizontal 
curvature and vertical curvature filter arrays with the following coefficient values: 
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Filter Array I: Horizontal Curvature Filter Array 
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Filter Array II: Vertical Curvature Filter Array 

By applying filter array I and filter array II to pixel array I [above), the 
following horizontal and vertical curvature estimates may be obtained: 

Horizontal Curvature - tt(A-4-B + C + 2"D + 2E + F-4-G + H) (14) 
Vertical Curvature - VWA + 2 B + C - 4*D <-4*E + F + 2-G + H) (IS) 

Thus, in this embodiment, curvature values are estimated based upon the values 
of pixels surrounding the deiect pixel, 

Still other embodiments are within the scope of the claims. 
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WHAT IS CLAIMED IS: 

1. A method of correcting one or more defect pixels (X) in a. source ^ 
image formed from a plurality of pixels each having a respective value, 
comprising correcting a defect pixel (X] based upon curvature information 
computed from pixel values located near the defect pixel (X). 

5 2. The method of claim 1, wherein the defect pixel (X) is corrected 

based upon curvature information computed for two different image directions. 

3. The method of claim 1, further comprising computing a curvature 
value based upon values of pixels surrounding the defect pixel (X). 

4. The method of claim 1, wherein the defect pixel (X) is corrected by 
10 interpolating values of pixels aligned with the defect pixel (X) in an interpolation 

direction determined from the curvature information. 

5. The method of claim 4, wherein the defect pixel (X) is assigned a 
value computed from a linear interpolation of values of pixels aligned with the 
defect pixel (X) along the interpolation direction, 

is S, The method of claim 5> wherein the defect pixel (X) is assigned a 

value computed from a linear interpolation of values of two non-defective pixels 
located on opposite sides of the defect pixel (X] > 

7. The method of claim 6, wherein the assigned defect pixel value is 
computed from the operations of adding (86) and shifting (88, 90). 

20 8. The method of claim 1 , further comprising obtaining the source 

image from an image sensor (20) having a plurality of light detectors each 
generating a pixel value of a respective color, wherein the colors are arranged in 
multiple color planes (40, 42, 44, 46) of non-overlapping pixel locations, and the 
defect pixel (X) is corrected based upon curvature information computed from 

25 values of pixels in the same color plane (40, 42, 44, 46). 



9. A method of correcting one or more defect pixels (X) in 'a source 
image formed from a plurality of pixels each having a respective value and 
corresponding to a respective color plane (40, 42, 44; 46), comprising: 

determining a median pixel value (M) from values of pixels on the same 
color plane (40, 42, 44, 46) located near a defect pixel (X); and 

correcting the defect pixel (X) based upon the median pixel value [M)> 

10. The method of claim 9, further comprising computing a difference 
value representative of the difference between the defect pixel value and the 
computed median pixel value (M), wherein the defect pixel (X) is corrected by 
replacing the defect pixel value with the median value (M) in response to a 
determination that the difference value exceeds a threshold value (T), 
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